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ABSTRACT 

Generalized Scidar (GS) measurements taken at the Paranal Observatory in 
November/December 2007 in the context of a site qualification for the future European 
Extremely Large Telescope E-ELT are re-calibrated to overcome the bias induced on 
the C^f profiles by a not correct normalization of the autocorrelation of the scintilla- 
tion maps that has been recently identified in the GS technique. A complete analysis of 
the GS corrected measurements as well as of the corrected errors is performed statis- 
tically as well as on individual nights and for each time during all nights. The relative 
errors of the C\ profiles can reach up to 60% in some narrow temporal windows and 
some vertical slabs, the total seeing up to 12% and the total integrated turbulence 
J up to 21%. However, the statistic analysis tells us that the absolute errors of the 
median values of the total seeing is 0.06 arcsec (relative error 5.6%), for the boundary 
seeing 0.05 arcsec (relative error 5.6%) and for the seeing in the free atmosphere 0.04 
arcsec (relative error 9%). We find that, in spite of the fact that the relative error 
increases with the height, the boundary and the free atmosphere seeing contribute 
in an equivalent way to the error on the total seeing in absolute terms. Besides, we 
find that there are no correlations between the relative errors and the value of the 
correspondent seeing. The absolute error of the median value of the isoplanatic angle 
is 0.13 arcsec (relative error 6.9%). 
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1 INTRODUCTION 

In November/December 2007 a site testing campaign (here- 
after PAR2007) has been organized at Cerro Paranal by ESO 
with the collaboration of several international teams (Dali 
Ali et aL]|2010 |. A set of different instruments run simulta- 
neously for a number of nights variable between 10 and 20 
nights with the aim to quantify different parameters charac- 
terizing the optical turbulence (the seeing e, the isoplanatic 
angle 6q, the wavefront coherence time To, the spatial coher- 
ence outer scale Co , the vertical profile of the spatial coher- 
ence outer scale Co(h), ....). The instruments employed were: 
a DIMM (Differential Image Motion Monitor), a GSM (Gen- 
eralized Seeing Monitor), a MASS (Multi Aperture Scin- 
tillation Sensor), a LuSci (Lunar Scintillometer), a MOSP 
(Monitor of Outer Scale Profiler) and a Generalized Sci- 
dar (GS). Some of these instruments (the DIMM and GSM) 
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provide integrated values of the turbulence energy (e, 9o, 
to, etc.) developed above the ground up to 20-25 km. Some 
of the instruments (the GS, the MASS, LuSci) are vertical 
profilers of the turbulent energy {C% profiles) or other astro- 
climatic parameter such as the spatial coherence outer scale 
(MOSP). One of the main goals of the PAR2007 site testing 
campaign was the inter-comparison of measurements that 
is a crucial passage not only for the characterization of a 
site, but, even more, for a characterization of the accuracy, 
uncertainty and dispersion of each instrument with respect 
to others. Dali Ali et al. (20101 presented a summary of 



the analysis of all the measurements retrieved from all these 
instruments from 17 December up to 26 December 2007. 

The GS is a vertical profiler based on an optical remote 
sensing technique that is able to reconstruct the turbulence 
distribution developed on the whole atmosphere from the 



ground up to 20-25 km. The vertical resolution (Vernin & 



Azouit (1983 1 ) AH increases with the height h and it de- 
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Table 1. Binary stars observed during the PAR2007 site testing campaign (20 nights) extracted from the Bright Stars Catalogue. In 
the first column the name of the star, in the second and third columns the right ascension and the declination, in the fourth column the 
magnitude of the primary star, in the fifth column the difference in magnitude between the primary and secondary star of each binary, 
in the last column the angular separation. 



pends on the binary separation, the wavelength and the ge- 
ometrical optical set-up as: 



AH(h) 



0.78 • y/\\h-h gs \ 



(1) 



where h is the height from the ground, h gs is the con- 
jugated height under-ground (see Section [5| and 6 is the 
binary separation. Considering the typical values of the bi- 
naries used for the GS, we can say that the GS vertical 
resolution is typically of the order of 1 km. The GS is, at 
present, one of the most useful and reliable instrument for a 
deep and detailed characterization of the optical turbulence 
and its spatial distribution in the atmosphere because, with 
the C% profiles and the wind speecQprofiles, it is possible to 
retrieve basically all the key parameters characterizing the 
optical turbulence for astronomical applications. The prin- 
ciple of the GS is solid, simple and many GS instruments 
have been developed in the last decade by different teams 
in the world. However, the GS is not a practical instrument 
for long terms and systematic monitoring of the turbulence 
and its distribution above an astronomical site because it re- 
quires a telescope with a pupil size of at least 1 m. For this 
application it is preferable to use automatic monitors such 
as the MASS, for example. The latter has a lower resolution 
(A/i ~ h/2) even if, recently, solutions to increase the reso- 
lution in the troposphere have been proposed - ( |Kornilov fc| 
Kornilov|20l"Tj ). 

The necessity to recalibrate the GS measurements of the 
PAR2007 campaign appeared urgent after the recent stud- 
ies \ Johnston et al.|(|2002|| and | Avila et al.|(|2009|) that put 



in evidence an error in the procedure for retrieving the C N 
profiles in the standard GS technique (in the post-processing 
phase) that could induce an overestimate of the turbulence. 
The error is perfectly described analytically therefore mea- 
surements can be corrected in a post-processing phase. In 



1 The wind speed can in principle be retrieved from the GS that 
has been applied with this goal in many studi es ||Avila et al. 



(|2001|),|Garcia-Lorenzo et al.|(|2006),| Avila et al . (2006), Egneret 



|al.| ( 12007} , |Garcia-Lorenzo et al.| p009\ , |Masciadri et al.|(|2010}) 

However, it is worth noting that it has been proved l |Masciadri| 
|fc Garfias| l |2001| , |Hagelin et~al] ( |2010[ )) that, wind speed pro- 
files from meso-scale models can supply equivalent reliable esti- 
mates of this parameter with the advantage of a larger spatial 
and temporal cover. For a systematic turbulence monitoring this 
is a preferable solution. 



this paper we summarize the results obtained after the re- 
calibration of the data-set. A statistical analysis as well as 
an analysis on each single night and for each instant along 
all the nights is performed so to have an absolute unbiased 
reference for the GS measurements. 

In Section [2] we describe the main principle of the GS 
technique, where the error of the GS technique has been 
identified and how to correct it. Effects of the errors in- 
duced on the vertical C% profiles are discussed in Section 
|3.1| Effects of errors induced on the total seeing as well as 
the boundary and free atmosphere contributions are treated 
in Section 13.21 In Section 13.31 we treat the effects of errors 
induced on the isoplanatic angle. Besides, in Section [3] a 
statistic of the median values of total seeing, boundary and 
free atmosphere seeing as well as isoplanatic angle obtained 
after re-calibration are provided. A brief discussion is pre- 
sented in Section[4] In Section[5]the conclusions of this study 
are presented. 



2 THE INSTRUMENT AND THE 
RE-CALIBRATION 

The SCIDAR technique (Scintillation Detection and Rang- 



ing) has been originally proposed by Vernin & Azouit ( 1983 1 



and Rocca et al. ( 1974 I to measure the vertical optical tur- 



bulence distribution (the refractive index structure constant 
C% profiles) in the troposphere. The technique relies on the 
analysis of the scintillation images generated by binary stars 
on the pupil plane of a telescope. The standard SCIDAR 
technique (called Classic Scidar) is insensitive to the tur- 
bulence near the ground. This fact represented in the past 
an important limitation for monitoring turbulence for as- 
tronomical applications because it is known that most of 
the turbulence develops in the low part of the atmosphere. 
To overcome this limitation ( jFuchs et al| ( |1994[ ), |Fuchs et] 
al. (1998)) proposed a generalized version of the SCIDAR 
(called Generalized SCIDAR - GS) in which the detector is 
virtually conjugated below the ground at a distance h gs per- 
mitting to extend the measurements range to the whole at- 
mosphere (from the ground up to ~ 20-25 km). GS measure- 
ments have been later done by several authors above differ- 
ent astronomical sites using instruments developed by sev- 



eral different teams (Avila et al. (19971, Avila et al. (19981, 



Avila et al. (20041, Kluckers et al. 


(19981, Mc Kenna et al. 


(20031, Garcia-Lorenzo et al. (2006 


I.Egneretal. (20071, Eg- 



© 2011 RAS, MNRAS 000,[T|fTol 



Re- calibrated Generalized- Scidar measurements at Cerro Parana! 3 





Figure 1. Left: pupil of the Auxiliary Telescopes (ATs) of the VLT at Cerro Paranal: £1=1.8 m, e=0.07. Right: autocorrelation of the 
telescope pupil shown on the left panel. 



et al 



et al. 



ner & Masciadri (20071, Fuensalida et al. (20081, Masciadri 



20101, Garcia- Lorenzo et al. (20111, Garcia-Lorenzo 



2011b I, Avila et al. (2011 1) 



The GS PAR2007 measurements were taken from an 
Auxiliary Telescope (AT) of the VLTI with CUTE-Scidar 
III instrument, which has been developed by the Instituto 
de Astrofsica de Canarias (IAC) team (Vazquez-Ramio et 
al., 2008) to be used at the Paranal Observatory within the 
European project ELT design study. This instrument is an 
improved version of others previously developed by the same 
group for Roque de los Muchachos and Teide Observato- 
ries, which has been used since 2002 to extensively monitor 
the optical turbulence at both observatories (Fuensalida et 
al. 



2004, 



Garcia-Lorenzo et al. (20111, Garcia-Lorenzo et 



al.| ( |2011b| )). CUTE-Scidar HI instrument provides pro- 
files in real time without dome and mirror seeing compo- 
nents. The binary stars specifications observed during the 
PAR2007 site testing campaign are listed in Table [T] Con- 
sidering the combinations of values of 6 and h gs used in 
the observations, we calculated that, using Eq|3] the typical 
vertical resolution AH at h=0 is within 400-1000 m. The 
AH(h) increases as Eq|3] 

The GS is based on the observation of binaries having 
an angular separation typically of ~ 3-10 arcsec, magnitude 
< 5-6 mag and a Am ~ 1 mag. When two plane wavefronts 
coming from a binary and propagating through the atmo- 
sphere meet a turbulent layer located at a height h from the 
ground, they produce, on the detector plan, optically placed 
below the ground at a few kilometers (h gs ), two scintillation 
maps characterized by typical shadows appearing in couple 
at a distance d. Such a distance is geometrically related to 
the position of the turbulent layer as d = 9(h+h gs ). The 
calculation of the auto-covariance (AC) of the scintillation 
map, normalized by the average of the intensity of the scin- 
tillation maps, produce the so called 'triplet'. The central 
peak is located in the centre of the AC frame; the lateral 
peaks are located at a symmetric distance d from the cen- 
tre. The amplitude of the lateral peaks is proportional to the 
strength of the turbulence of the layer located at the height 
h weighted by the scintillation that such a layer produces 
on the detector. In a multi-layers atmosphere, different tur- 
bulent layers (i) produce triplets with lateral peaks located 
at different distances di from the centre of the AC frame. 



The C% profiles are obtained inverting the normalized func- 
tion just described called Friedholm equation that describes 
several triplets. 

In 2002, Johnston et al. put in evidence an error in the 
calculation of the normalization of the AC frames and they 
studied the effects of this error on the C% profiles for h 
= 0. More recently, Avila & Cuevas (2009), extended the 
analysis for heights h > 0. Results of this analysis say that, 
to obtain exact results for the C% profiles at all heights h 
from the ground, one has to multiply the C%(h) retrieved 
from the GS by a factor (l/l+£(ft)) where £ is the relative 
error between the exact and erroneous auto-covariance of 
the scintillation maps obtained with the GS and described 
in Eq.|2}: 



C(A) 



S(r-eh gs ) 



aS(r) + b[S(r + 6h gs ) + S{r - 6h g3 )] 



(2) 



where S(r) is the auto-correlation of the pupil of the tele- 
scope, a and b are related to Am as Eq.(|3|: 



1 + 
+ 



b = 



(l + o) 



a = 10" 



(3) 



C(h) depends on a set of geometrical parameters related 
to the optical set up and the observed binaries, more pre- 
cisely the pupil of the telescope D, the height h gs at which 
the detection plane is conjugated below the ground, the ratio 
between the stellar magnitudes of the binaries Am, the an- 
gular separation of the binary 9 and the ratio e between the 
central obscuration D* and the telescope pupil (e=D* / D). 
We obtain therefore £=£(0,h gs ,e,b,h). Figure [I] shows the 
pupil size of the Auxiliary Telescopes of the Very Large Tele- 
scope where the GS has been run during the PAR2007 cam- 
paign and a section of the bi-dimensional autocorrelation of 
the pupil of the telescope S. We note that the value of e for 
the AT at Paranal is smaller than in other telescopes used 



for GS measurements (Masciadri et al. (20101 Avila et al 



( 2011 1 |Garcia-Lorenzo et al. (20111). As it will be discussed 



in Sec HI this is a favorable condition in the context of the 
error of the normalization of the scintillation maps of the 
GS. 
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Figure 2. Relative errors of the vertical profiles calculated for all the 20 nights monitored during the PAR2007 site testing 

campaign. The titles of the panels report the name of the observed binary stars with the correspondent angular separation. In each 
panel the relative errors with respect to the height h associated to different values of the h gs (i.e. the height at which the CCD plane is 
conjugated underground) are shown. The numbers in parenthesis indicate the percentage of profiles with respect to the whole sample of 
measurements . 



3 RESULTS 

3.1 Vertical turbulence distribution: C% 

Figure [2] reports the relative errors ((h) of the C% profiles 
(before and after re-calibration) as a function of the height h 
calculated for the 20 nights monitored during the PAR2007 
site testing campaign. For each binary star (i.e. for a fixed 
angular separation 9 and magnitude ratio Am) the rela- 
tive error is calculated for different values of the conjugated 
height under the ground (h gs ) (full, dashed and pointed 
lines). Panels in Fig|2] include all the combinations of the 
set of parameters from which the error ( depends on. Typi- 
cally 3 or 4 binary stars are observed in each night and the 
position of the conjugated plane has been changed sometime 
during the nights. As a consequence, in each night we have 
the sequence of a few of the ((h) shown in Fig[2] We observe 
in Fig(2]that ( assumes a value of the order of 10% for h 



=0 in almost all of the cases but it increases monotonically 
with the height reaching also values of the order of 40-60% 
in some cases above 15 km (BS1563). It is however worth 
noting that this happens for a small percentage of cases (see 
numbers in parenthesis in the panels of Figj2|. Indeed, if 
we look at the average C% profiles obtained in each night 
before and after the re-calibration (Annex [A|, we note that 
the vertical distribution (shape of the C%) remains substan- 
tially the same for each height and the difference between 
the original and the corrected C% profiles appears small. 
This fact indicates that it is important to consider always 
both elements to appreciate the impact of the ((h) error on 
GS measurements. Looking at FigjX] we note that the shape 
of the C% profile is always characterized by a principal peak 
near the ground with most of the turbulent contribution. 
However, in the free atmosphere, we observe very different 
spatial distributions of the turbulence in each night and it 
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Figure 4. Relative error versus total seeing (left), seeing below (centre) and above (right) 1000 m for all the nights of the PAR2007 site 
testing campaign. Y-axis values have to be multiplied by 100 and expressed in (%). Ex: 0.1 means 10 %. 




0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 0.2 0.4 0.6 0.8 1 

Seeing Tot. (arcsec) Seeing h<1km (arcsec) Seeing h>1km (arcsec) 

Figure 5. Original (x-axis) versus re-calibrated (y-axis) seeing values for the 20 nights monitored during the site testing campaign of 
the PAR2007 site testing campaign. Each point is the average value of a single night. Left: total seeing. Centre: boundary layer seeing 
(h < 1km). Right: free atmosphere seeing (h > 1km). The error bars represent the standard deviation. 




Figure 3. Median vertical profile of the calculated on the 
whole re-calibrated data-set of the 20 nights of the PAR2007 site 
testing campaign. Dashed lines represent the first and third quar- 
tiles. 

is not necessarily visible the typical secondary peak at the 
jet-stream level that appears when we estimate the average 
or the median values of a rich statistical sample. This tells 
us that in the free atmosphere there is a large variability of 
the position of the layers in each night. As expected, in the 
re-calibrated C% profiles, the turbulence is systematically 
weaker than in the original ones at all heights. In Annex 
[B] are reported the temporal evolution of the re-calibrated 
C N profiles in each night of the PAR2007 campaign. Empty 



temporal windows correspond to an interruption of acquisi- 
tion on the same binary star in some cases or to the change 
of the binary star in some other cases. Looking at FigjBjand 
Fig|C] (Section [3T2] ) it is possible to identify the instants in 
which the binaries have been changed and it is possible to 
discriminate between the two cases. In case the binary star 
is changed, this means that the parameters from which the 
relative error (,(h) depends on are changed and therefore the 
value of the relative error changes sharply its value. We note 
that, most frequently, each night is characterized by a few 
layers (a sort of background structure) that are present dur- 
ing the whole night. Sometime it happens however, that a 
particular layer appears or disappears at a precise time dur- 
ing the night, for example on 10/11/2007 at 12 km at around 
04:15 UT; on 17/12/2007 at 13 km at around 06:00 UT; on 
19/12/2007 at 4 km at around 05:00 UT; on 20/12/2007 at 
2 km at around 06:00 UT. Depending on the height and the 
intensity of each turbulent layers these events can induce 
more or less rapid and important changes of the integrated 
astro-climatic parameters. For example, the disappearing of 
a layer at 15 km can frequently produce, a weak decreas- 
ing of the total seeing (at this height the turbulent layers 
are in general much weaker than near the ground) but an 
important increase of the isoplanatic angle that is particu- 
larly sensitive to even small turbulence changes in the high 
atmosphere. 

Figure [3] shows the median C% profile obtained with 
the re-calibrated data-set associated to all the 20 nights of 
the PAR2007 campaign. In the free atmosphere we identify 
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three major secondary peaks at 4 km, 7- 
14-15 km. 



km, 10-11 km and 



3.2 Integrated turbulence: seeing and J 

In Annex [C] are reported the temporal evolution of the rel- 
ative errors of the total seeing Etot, the boundary ebl 
and the free atmosphere seeing efa (before and after re- 
calibration) for all the 20 nights of the PAR2007 campaign. 
The relative error for the total seeing is defined as: 



ttot 



e TOT ~ e TOT 



-TOT 



where e tot and etot are: 

' 20km 

= -1(A)- i J C%{h) 





3/5 



dh 



20km 



3/5 



£tot — 



C%(h)dh 



(4) 



(•») 



(6) 



and A(\) is equal to 19.96xl0 6 for A=0.5 /xm. Similar equa- 
tions have been built for the relative error for ebl and efa 
replacing the extremes of the integrals with the ranges: [0- 
1000m] and [1000m,20km]. Looking at Figjc] we see that, in 
all cases, the relative error of the boundary layer is relatively 
modest, most of time of the order of 5-9%. The relative error 
of the free atmosphere seeing is slightly greater but always 
smaller than 15%. It is greater than 10% only in 5 occasions 
(16/11/2007, 17/11/2007, 14/12/2007, 16/12/2007) and for 
a partial part of the night. We note that, for the bound- 
ary layer contribution, the relative error is almost constant 
and it changes its values basically when the binary star is 
changed. On the other side, the relative errors of the free 
atmosphere fluctuate in time (sometimes in a not negligi- 
ble way but always more than in the boundary layer where 
the relative error remains basically constant in each seg- 
ment associated to each observed binary). Why ? It is true 
that the absolute error (Eq[2| increases with the height but 
this means that the relative errors at each instant and the 
successive one are both larger in the high part of the atmo- 
sphere than in the low part of the atmosphere. Therefore 
the fact that we note a systematic and more evident vari- 
ation of the relative error in successive instants in the free 
atmosphere and not in the boundary layer can not be re- 
lated to this reason. This is highly probably due to the fact 
that the C% in the free atmosphere is, in absolute terms, a 
couple of order of magnitude smaller than the C% near the 
ground. As a consequence, when we calculate the relative 
error in the free atmosphere, we find that small variations 
in the C% are associated to larger difference between the 
relative error at an instant to and the successive one. In 
other words, the variation of the numerator of Eqj4] at two 
successive instants is much more smaller than the denom- 
inator in the case of the boundary layer than in the case 
of the free atmosphere. Figure [4] shows the relative errors 
for the total seeing, the seeing in the boundary layer and 
in the free atmosphere as a function of the corresponding 
values of seeing for all the nights of the PAR2007 site test- 
ing campaign. It is evident that no correlation exist between 




1.5 2 
Seeing (arcsec) 

Figure 6. Cumulative distribution of the total seeing (black line), 
seeing below (red line) and above (green line) 1000 m calculated 
before and after the re-calibration. The re-calibrated seeing (thin 
style line) is correctly always smaller than the original one (bold 
style line). The numbers in the panels are the median values cal- 
culated before (bold style line) and after (thin style line) the 
re-calibration. 



the relative error and the corresponding value of the seeing. 
This means that it is not possible to retrieve any relation 
between the value of the relative error of the seeing and 
the value of the seeing. However, it is possible to note that 
the fluctuations of the relative errors in the free atmosphere 
are larger than in the boundary layer and the distribution 
of the relative errors of the seeing assumes a sort of 'cloud 
shape'. Figure [5] shows the re-calibrated averaged values of 
the seeing calculated in each nights versus the original av- 
eraged values. As expected, a systematic decreasing of the 
estimated turbulence in the re-calibrated data is evident. At 
the same time, we observe that the difference on the mean 
values (before and after the re-calibration) is very small if 
compared to the standard deviation of each estimate. The 
greater is the seeing, the greater is the absolute difference 
between the average seeing before and after re-calibration. 
Figure [6] shows the cumulative distribution of the total see- 
ing, seeing in the boundary layer and in the free atmosphere 
before and after the re-calibration. The corresponding me- 
dian values reported in Figj6]tell us that the absolute errors 
for the median values of the seeing are: AesL=0.05 arc- 
sec, Aei?A=0.04 arcsec, Aetot~0-06 arcsec. In Annex |E| is 
shown the temporal evolution of the seeing (total, boundary 
layer and free atmosphere for each night. 

It is worth noting (Fig[7jl that the relative error of the 
total energy J (measured in m 1 ^ 3 ) defined as: 



C 2 N (h)dh 



is higher than that of the seeing and it can reach 20% 
in the free atmosphere, as well as in the total atmosphere. 
We distinguished between J and seeing because, in some 
applications, it can be preferable to treat the total energy 
instead of the seeing. J is indeed linear with respect to ad- 
dition while seeing is not. 
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Figure 7. Relative error versus total J (left), J below (centre) and above (right) 1000 m for all the nights of the PAR2007 site testing 
campaign. Y-axis values have to be multiplied by 100 and expressed in (%). Ex: 0.1 means 10 %. 




Figure 8. Relative error versus isoplanatic angle (left) and cumulative distribution of the isoplanatic angle (8o) before (bold style line) 
and after (thin style line) the re-calibration of the data (right). 



3.3 Isoplanatic angle 

In Annex[D]is reported the temporal evolution of the relative 
error of the iso-planatic angle (do) for all the 20 nights of 
the PAR2007 campaign. 8q is defined as: 



erature (|Avila et al.| (|2004|), |Masciadri et aL] (|2004[) , |Egner 



8 = 0.057- 




-3/5 



(h)dh 



(8) 



where h up corresponds to the top of the atmosphere (more 
frequently to the end of vertical profile). We observe that the 
features of the temporal evolution of the 8o relative error are 
similar to those of the free atmosphere seeing. This is not 
surprising because the isoplanatic angle is particularly sen- 
sitive to the turbulence in the high part of the atmosphere 
due to the presence of the factor h 5 ^ 3 in the integral. The 
9o relative error is never larger than 18% with the exception 
of a few minutes on 17/11/2007 in which the error was of 
the order of 18%. Figure Ill-left shows the 8q relative error 
versus the corresponding value of 8o related to all the 20 
nights of the PAR2007 site testing campaign. Also in this 
well as for the seeing, no correlation between the 
two quantities is observed. Figure fright shows the cumu- 
lative distribution of the iso-planatic angle before and after 
re-calibration. As expected, the median value of the distri- 
bution increases after the re-calibration (from 1.74 arcsec 
to 1.87 arcsec) because the total turbulence decreases. We 
highlight that 8o in Fig|8] is calculated using h up —20 km 
(i.e. 20 km from the ground). Previous studies in the lit- 



et al. (20071, Masciadri et al. (20101, Hagelin et al. (I2011I), 



Garcia-Lorenzo et al. (20111, Garcia-Lorenzo et al. ( 2011b[ |) 
performed the calculation of #o with a value of h up between 
20 and 22 km from the ground. To estimate the effect of the 
value of h up on 8o and to evaluate potential differences due 
this element, we calculated the integral using h up =22 km. 
We find, respectively, 1.73 arcsec and 1.87 arcsec. We con- 
clude, therefore, that the median value of 8o remains the 
same (1.87 arcsec) within one hundredth of arcsec in these 
data-set. In Annex[F]is shown the temporal evolution of the 
isoplanatic angle for each night. 



4 DISCUSSION 

From a qualitative point of view we can say that, as shown 
in Avila & Cuevas (2009), the greater is the distance of the 
conjugated plane under ground (h ga ), the greater is C(ft); 
the greater is the binary separation 8, the greater is 
the greater is Am, the greater is C(/i); the smaller is the 
pupil size D, the greater is C(/l); the greater is the ratio be- 
tween the central obscuration and the pupil size e=D* /D, 
the greater is near the ground. What about the quanti- 
tative effects ? We put in the context our results with respect 
to others recently published. 

Masciadri et al. (20101, treating a sample of measure- 



ments related to 43 nights done above Mt. Graham (VATT 
telescope - D=1.83 m), proved that GS measurements ob- 
tained with a pupil size D > 1.83 m and a binary separation 
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< 8 arcsec are affected by this error for less than a few 
hundredths of an arcsec (0.04 arcsec corresponding to a rel- 
ative error of ~ 5% for the total seeing) with an average 
C% profile reproducing the same spatial distribution of the 
un-corrected profile. 

Similar results have been obtained by Avila et al. (2011) 
treating measurements done above the San Pedro Martir 
Observatory (D —1.5 m and D=2.1 m) in which it has been 
estimated a typical relative error of 4.2% (an absolute error 
of ~ 0.04 arcsec on the total seeing). Modest errors therefore 
even if, in that case, a binary with 8 = 14.4 arcsec (£ Uma) 
was selected among the targets. The modest relative error 
is due to the fact that the percentage of time of observation 
of £ Uma was relatively low (14.4 %). 

Results of this paper for measurements done at Cerro 
Paranal with a pupil size £>=1.8 m are coherent and similar 
to the previous ones. A relative error for the total seeing of 
5.6% (absolute error equal to 0.06 arcsec) and a relative error 
of the isoplanatic angle of 6.9% (absolute error equal to 0.13 
arcsec) are calculated. In just one case a binary separation 
larger than 10 arcsec has been used (BS1563 - 0=12.3 arcsec 
- percentage of observation is 7%). 

Some larger errors have been observed in the study per- 
formed at Roque de los Muchachos Observatories on Canary 



Islands ( Garcia-Lorenzo et al. (2011b l) because a set of dif- 
ferent geometric parameters (wider binary separation and 
deeper conjugation heights of the detector plane) have been 
used in combination with a smaller telescope pupil size. In 
the study of Roque de los Muchachos (D—l m), five binaries 
with angular separation wider than 10 arcsec (for a percent- 
age of time equal to 3.39%) have been selected. The absolute 
error of the median values of the total seeing is of the order 
of 0.12-0.15 arcsec and that of 9 of the order of 0.34-0.92 
arcsec (for observations at high and low vertical resolution). 
We conclude that, considering the relative small percentage 
of observations with wide binaries, the element that mostly 
affected the overestimate of the turbulence in the study done 
at Roque de los Muchachos was the small size of the pupil 
of the telescope (D=l m) and a deep conjugation heights of 



the detector plane. However, looking at Fig. 5 (Avila et al 



12009)), we observe that, from a quantitive point of view, 



the impact of the telescope pupil size is much more impor- 
tant than the conjugation heights of the detector plane. The 
size of the telescope seems therefore to be the most critical 
parameter in the error of the normalization of the average 
of the scintillation map in the GS measurements published 
so far. To contain the relative error it is suggested to keep 
the value of the conjugated height under ground h gs at 2 or 
3 kilometers maximum and to use a telescope with a pupil 
size of at least 1.5 m. The binary separation affected less 
the statistical results published so far because the number 
of measurements obtained with wide binaries was relatively 
small. 



tion of the optical turbulence (C% profiles) indicates relative 
errors that can achieve order of 40-60% (Fig[2| in some part 
of the atmosphere and in some narrow temporal windows. 
However, when we calculate the relative errors of the in- 
tegrated counterparts, the relative error reach a maximum 
value of 13% for the correspondent total seeing Etot and 
15% for the free atmosphere. The absolute errors are much 
smaller, particularly if one considers the C% average on a 
whole night. The shape of the average of the C% profiles 
before and after re-calibration is basically the same. For 
what concerns the integrated astroclimatic parameters, from 
a statistical point of view, we find a relative error of the me- 
dian value of the total seeing equal to 5.6%, for the seeing 
in the boundary layer equal to 5.6% and for the seeing in 
the free atmosphere seeing equal to 9%. These values cor- 
respond to absolute errors respectively equal to 0.06 arcsec, 
0.05 arcsec and 0.04 arcsec. The absolute errors provided 
by the boundary layer and the free atmosphere contribu- 
tions are statistically comparable even if the relative error 
increases with the height. The relative error of the turbu- 
lent energy J is slightly larger. It can reach a maximum 
of 21% in the total atmosphere and 23% in the free atmo- 
sphere. The relative errors for the isoplanatic angle are never 
larger than 18%. The absolute error of the median value of 
9o is 0.13 arcsec. The re-calibration of the GS data-set per- 
formed in this paper suggests a revision of results obtained 
in ( |Dali Ali et al.|[2010[ > on the cross-comparison between 
optical turbulence measurements obtained with different in- 
struments during the PAR2007 campaign. A few of us are, 
at present, deeply involved in a cross-checking of such a mea- 
surements with predictions performed with an atmospheri- 
cal non-hydrostatic numerical model and a forthcoming pa- 
per will be dedicated to this topic. 

Results of this paper, joint with results found on pre- 
vious papers on the same topic, indicate that, if a GS is 
used with a pupil size D > 1.8 m and an angular separation 
smaller than 10 arcsec (preferably with a value of the con- 
jugated plane not deeper than ~ 3 km underground), the 
absolute errors can be considered negligible (of the order of 
a few hundredth of arcseconds) from a statistical point of 
view. 

If D < 1.8 m and/or the angular separation is greater 
than 10 arcsec, the re-calibration is highly suggested, in par- 
ticular if one is interested on some specific temporal win- 



dows. We remember that it has been proved (Masciadri et 
al"1|2010| ) that the GS-HVR technique |Egner fc Masciadri 
(20071 to be applied in first kilometer from the ground is 



not sensible to the error of the normalization of the auto- 
correlation. This solution can therefore be selected to in- 
vestigate the turbulence at high vertical resolution near the 
ground (first kilometer) without necessity to re-calibrate the 
measurements. 



5 CONCLUSION 

In this paper a re-calibration of the whole data-set of mea- 
surements performed with the Generalized Scidar (CUTE- 
Scidar) at Cerro Paranal on November/December 2007 
(PAR2007 site testing campaign) for 20 nights is performed. 
The statistical analysis of the spatial and temporal distribu- 
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APPENDIX A: 



AVERAGE C% 



PROFILES 



Fig Al and Fig 



A2 



show the average C N profiles for each 
night of the PAR2007 campaign before and after the data- 
set re-calibration. The dates are in local time (LT). 



APPENDIX B: 
EVOLUTION 



C 2 N PROFILES TEMPORAL 



Fig Bl and Fig 



B2 



show the temporal evolution of the C N 
profile during the night for all the 20 nights of the PAR2007 
site testing campaign. 



APPENDIX C: TEMPORAL EVOLUTION OF 
THE RELATIVE ERROR OF THE SEEING FOR 
THE INDIVIDUAL NIGHTS 

From Fig |Cl| up to Fig |C4| are shown the temporal evolution 
of the relative error of the total seeing, the boundary layer 
[h < 1000 m) and the free atmosphere (h > 1000 m) con- 
tributions for all the 20 nights of the PAR2007 site testing 
campaign. 



APPENDIX D: TEMPORAL EVOLUTION OF 
THE RELATIVE ERROR OF THE 
ISOPLANATIC ANGLE FOR THE INDIVIDUAL 
NIGHTS 

Fig |D1 1 shows the temporal evolution of the relative error of 
the isoplanatic angle for all the 20 nights of the PAR2007 
site testing campaign. 



APPENDIX E: TEMPORAL EVOLUTION OF 
THE SEEING FOR THE INDIVIDUAL NIGHTS 

Fig |El| and Fig |E2| show the temporal evolution of the seeing 
(total, boundary layer and free atmosphere) for all the 20 
nights of the PAR2007 site testing campaign. 



APPENDIX F: TEMPORAL EVOLUTION OF 
THE ISOPLANATIC ANGLE FOR THE 
INDIVIDUAL NIGHTS 

Fig |Fl| and Fig |F2| show the temporal evolution of the isopla- 
natic angle for all the 20 nights of the PAR2007 site testing 
campaign. 
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Figure Al. Average vertical profiles for each night of the PAR2007 site testing campaign obtained with the original data-set 
(full line) and the re-calibrated data-set (dashed line). Ex: 9/11/2007 corresponds to 10/11/2007 in UT. Error values are expressed in 
percentage. 
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Figure A2. It follows Fig.jAl] 
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Figure Bl. Temporal evolutions of re-calibrated vertical profiles for all the 20 nights monitored during the PAR2007 site testing 
campaign (units are m~ 2 / 3 and in log scale). The dates are in local time (LT). Ex: 9/ff/2007 in LT corresponds to 10/lf/2007 in UT. 
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Figure CI. Temporal evolution of the relative error of the total seeing etot (left), the boundary layer seeing ebl (centre) and the 
seeing in the free atmosphere Efa (right) during the 20 nights. Ebl is defined as the seeing integrated for h < 1km, Epa as the seeing 
for h > 1km. The dates are in local time (IT). Ex: 9/11/2007 corresponds to 10/11/2007 in UT. Y-axis values have to be multiplied by 
100 and expressed in (%). Ex: 0.1 means 10 %. 



© 2011 RAS, MNRAS 000,[T}{To] 



16 E. Masciadri et al. 



4>. 



10 12 



10 12 



10 12 



10 12 2 



10 12 



10 12 



20-1 1-2007 



20-1 1 -2007 



10 12 2 



10 12 



10 12 



21-1 1-2007 



21-1 1 -2007 



10 12 



10 12 



10 12 



23-1 1-2007 






23-1 1 -2007 




23-11-2007 






E ° 




E 





10 12 



10 12 



10 12 



Figure C2. It follows Fig. Kg] 
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Figure C3. It follows Fig.[CT1 
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Figure C4. It follows Fig. Kg] 
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Figure Dl. Temporal evolution of the relative error of the isoplanatic angle 9q during the 20 nights of the site testing campaign of 
November/December 2007. The dates are in local time (IT). Ex: 9/11/2007 corresponds to 10/11/2007 in UT. Y-axis values have to be 
multiplied by 100 and expressed in (%). Ex: 0.1 means 10 %. 
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Figure D2. It follows Fig.|DT] 
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Figure El. Temporal evolution of the seeing during the 20 nights of the site testing campaign of November/December 2007. The dates 
are in local time (LT). Ex: 9/11/2007 corresponds to 10/11/2007 in UT. Black line: total seeing. Green line: seeing in the boundary layer 
(h < 1000 m). Light blue line: seeing in the free atmosphere (h > 1000 m). 
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Figure E2. It follows Fig. [El] 
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Figure Fl. Temporal evolution of the isoplanatic angle during the 20 nights of the site testing campaign of November/December 2007. 
The dates are in local time (IT). Ex: 9/11/2007 corresponds to 10/11/2007 in UT. The integral is calculated up to 20km from the 
ground. See text for discussion. 
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Figure F2. It follows Fig. [FT] 
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